Spin-dependent Rabi oscillations in single quantum dot 
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Ultrafast optical pump of exciton resonance in single quantum dot by an elliptically polarized 
laser pulse is described by the non-Markov balance equations. Population and spin dynamics are 
investigated and spin-dependent Rabi oscillations are considered for the case of a flat quantum dot, 
with the lateral size greater than the height. The following peculiarities of temporal evolution have 
been found: a) a beating of Rabi oscillations under an elliptically polarized pump and b) quenching 
of spin orientation due to a spin-orbit splitting effect. 
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PACS numbers: 73.21. La, 78.47.+p, 73.20.-r 

Investigations of excitonic excitations in a single quan- 
tum dot (SQD) have been started in the past decade 
(see recent Refs. in |]J). The Rabi oscillations and spin- 
related effects in a QD ensemble have also been inves- 
tigated, see 0] and 0] respectively. Recently the Rabi 
oscillations of electron population has been reported in 
SQD excited by ultrashort laser pulse UQ- Due to this, 
an examination of the spin-dependent ultrafast response 
in SQD is now timely. In this paper, we consider the spin- 
dependent coherent dynamics of an exciton in a SQD 
under ultrafast interband pump by elliptically polarized 
radiation. 

A temporal evolution of the exciton is described by 
the 2x2 spin-dependent density matrices written in the 
basis of the vacuum and excited excitonic states. These 
density matrices are governed by the balance equations 
with a non-Markov generation rate which is dependent on 
the population redistribution. For the simplest case of a 
weak spin-orbit interaction, allowing the only interband 
transitions between the same spin states, a + or — circu- 
lar polarized pump can only excite spin-down or spin-up 
states. Since a temporal response under an elliptically 
polarized pump can be considered as a superposition of 
Rabi oscillations under ± polarized pumps with differ- 
ent intensities, a beating of Rabi oscillations takes place. 
For the more complicated case involving spin-orbit in- 
teraction in c-band, one can obtain a quenching of spin 
orientation in a SQD even under the circular polarized 
pump. Below, we have considered these phenomena in 
the framework of a flat QD model, with the lateral size 
greater than the height, and have estimated the maximal 
degree of spin orientation. 

We consider the coherent dynamics of a QD under res- 
onant interband excitation described by the perturbation 
operator [SH t exp(— iu)t) + H.c.]. The averaged over pe- 
riod, 27T /u>, density matrix, R t , obeys the following quan- 
tum kinetic equation (see evaluation in 



dRt 
dt 
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where 5 — ► +0 and S t = exp(iHt/~h) describes evolu- 
tion of QD with the Hamiltonian H. The transitions 
between the spin-degenerate c-band and v-band states 
are described by the matrix 



SHt 



-Ew t 



g e 
gt o 



(2) 



Here we consider the interband photogeneration caused 
by the electric field E[e exp(— iu>t) + c.c.]wt, with the unit 
complex polarization vector e, and Wt = exp[— (i/r p ) 2 ] 
is the Gaussian form-factor of the pulse, r p stands for 
the pulse duration. The polarization-dependent matrix 
g e = (e • v cv ) is written through the interband velocity 
matrix element, ~v cv , and g e is the non-Hcrmitian matrix 
because e ^ e*. 

Below, we shall restrict our consideration to the vicin- 
ity of the excitonic resonance when R t can be replaced 

Pe,t 

Po,t 

and p 0jt describing the excited (e-) and vacuum (o-) 



by the diagonal matrix 



with 2x2 matrices 



Pe,t 



states respectively. In agreement with the experimen- 
tal conditions , we have neglected here the two-photon 
biexcitonic resonance, see 0)13 an d Refs. therein. Thus, 
we transform Eqs.(l,2) into the following form: 



0_ 

dt 



pe,t 
Po,t 



TlbJ 



J — OO 

x {g e p ,t> - Pe,t>g e ) g~t + H.c, (3) 



where hAu is the detuning energy with respect to the 
exciton peak energy. Adding these Eqs. we obtain that 
Pet + Pot is a time-independent matrix which is deter- 
mined through initial conditions (see below). 

To calculate the matrix g e we neglect the Coulomb 
renormalization of the interband velocity matrix element 
and consider 



<A|v|A') 



33' 



jj' 



(4) 



Here v^y is the 8x8 interband velocity matrix and x is 
the in-plane vector. For the flat QD model, we write 



2 



the overlap integral in (3) through the in-plane enve- 
lope functions ipj X - Here we have replaced the transverse 
contribution to the overlap integral by unit considering 
transitions between top v-band states and bottom c-band 
states. Neglecting the nonparabolic (oc p 3 ) contributions, 

we describe the w-band states by (/j x \a), where tp^ h ^ 
takes into account the heavy- light hole mix [6( and \a) 
corresponds to the degenerate spin state directed along 
OZ with the spin numbers a = ±1. Taking into account 
noticeable spin-orbit effects in c-band we write the eigen- 
state problem for spinor ?/> x ^ as follows: 



Im 



[v s xp])- £ U«W=0, 



(5) 



where <r are the Pauli matrices, m is the electron effective 
mass, and the characteristic spin velocity v s is along the 
growth axis, OZ 0,0. 

We shall restrict our consideration to the case of weak 
spin-orbit effect, performing the unitary transformation 

— exp[— im(cr ■ [v s x x\)/h]ip x . Taking into account 
oc v 2 contributions only, we obtain the Schrodinger equa- 
tion in the form: 



(ff^xx])/ft]-eU = 0, (6) 
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We consider this equation with the zero boundary con- 
dition tpjfp = 0; the curve r is bounded by a flat QD 
with the area S. Note that the momentum operator here, 
[p x x\/h, is along OZ, so that the spinor ip x is propor- 
tional to the above introduced spinor \a). 

Next, using the two-level basis including the lower elec- 
tron states with a = ±1 and the upper hole states with 
a' = ±1, we obtain the only non-zero matrix elements (2) 
for the transitions \e, 1) — > Jo, 1) and \e, — 1) — * \o, — 1); 
see similar selection rules in For the case of ellipti- 

cally polarized excitation described by e = A + e + + A_e_, 
where e± are the polarization vectors for ± circular waves 
and A± are complex numbers (moreover |A+| 2 + |A_| 2 = 
1), we obtain the matrix g e in the form: 



VfiJ 



A_ 
-X A 



The non-zero coefficients here are written through 



v fi± = / dx^V x W 



1 



/ mv s \ - 



s 



(7) 



(8) 



where V is the Kane velocity. The overlap integral here 
have been calculated for the cases of circular and square 
QDs when the numerical coefficient a is equal to 0.11 and 
0.32 respectively. A suppression of the matrix element (7) 
with increasing v s also takes place for the conic QD case 

m 

We describe the coherent dynamics by the use of the 
diagonal population numbers for the excited and ground 
states, n e ,<ri = (e, <r\p e t\e, a) and n . at = (o,a\p 0<t \o,o). 



The system (3) is transformed into the independent equa- 
tions for the states with cr = ±1: 



d 
dt 



n e fi±it 
n ,&±it 
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Here we have introduced the kernel <&* 



(9) 



wtWt' cos Auj(t — t'), and v r 



2{eE/Tu J j) 2 \v 6A 



is the characteristic photogeneration rate. Since 
n e ,a±it + n ,B±it are independent of time, and the 
oscillation strengths for ± transitions are proportional 
to |A T | 2 , we use the phenomenological conditions in the 
form 



«e,6±lt + ™o,fi±lt — |A T | 2 , 



(10) 



which is in agreement with the particle conservation re- 
quirement z\^ a {n e at + n oa t) = 1. Using Eqs.(9), we ob- 
tain the closed equation for the population differences, 
An± t = n a e±it - n e «±it, in the form: 



dAn± t . |2 



dt' 



$ tt ,An ±t = (11) 



= | A=p | 2 which are 



with the initial conditions An±t_> 
obtained from Eq.(10) and Ti e ,B±it— >- <x> = 0. 

The concentration of photoexcited electrons is given 
by nt = Ti e fiu + n e fi-\t while the spin orientation, St, is 
determined as follows St = X) CT =±i( e ' ^I^K a ) n e,at- The 
only z-component of spin appears to be non-zero and 



k(n e ,m - riefi-it) = kSn t , 



( mv s \ ' 



s. 



(12) 



where k = J d~x.ip x ^ + a g ipx and a is the coefficient from 
Eq.(8). The spin orientation is proportional to the redis- 
tribution between the ± spin states, 5nt, and the factor 
k decreases with increasing v s . 

For the resonant excitation case, Aw = 0, the kernel 
in Eq.(ll) is not dependent on t — t', one obtains the 
second-order differential equation 

fAn^ _ 2t dAn^ + •V,,,^,/^ = „ 



dt 2 



dt 



.(13) 

with the use of an additional initial condition 
[w^ 1 dAns±t/dt]t~>-oo — 0. The solution of Eq.(13) takes 
the form: 



An B± t = \X T \ cos 



A\X, 



t/7 



dre 



(14) 



where A = ^Txv r T v is the excitation amplitude and 
|A±| 2 = (1 ± s)/2 with the degree of circular polariza- 
tion s = |A + | 2 
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FIG. 1: Temporal evolution of the total population, n t , and 
the spin redistribution, Sn t , under the excitation amplitudes, 
A — 7T (a) and A = 3.5tt (b). Solid, dashed and dotted curves 
correspond n t under circular, elliptical and linear polarized 
pumps respectively. The dot-dashed curve presents Snt under 
elliptical polarized pump. The pulse profile is shown as the 
short-dot curve. 



In Fig.l we plot the temporal evolution of n t and 5n t 
for the cases of circular, elliptical (with s = 1/2), and lin- 
ear polarized pumps under the resonant excitation con- 
dition Aiv — 0. Note that 5n t = —n t and Sn t — for the 
cases of circular and linear polarized pumps respectively. 
For the 7r-pulse excitation case (if A = tt), one obtains 
monotonic increase of nt with nt^oo = 1 for the circu- 
lar pump. The spin redistribution, 8n t < 0, varies with 
t/r p in a similar manner. As the excitation amplitude 
increases, the temporal evolution becomes oscillatory, as 
shown in Fig. lb for A = 3.57T. These dependencies have 
different periods, for the circular, linear, and elliptically 
polarized pumps. 

Next we consider the detuning effect by calculating 
the solution of Eq.(ll) numerically. In Fig. 2 we plot the 
temporal evolution of n t and Sn t under different Awt p . 
The detuning results in suppression of both the Rabi flop 
after the pulse and the amplitude of oscillations if A > n. 
The characteristics of the temporal response for different 
polarizations of pump do not change in comparison with 
the resonant case. 

The oscillations of Rabi flopping population, n = 
n t ^oo, and the redistribution between the ± spin states, 
Sn = 5n t ^oo, are shown in Fig. 3 for the resonant ex- 
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FIG. 2: Quencing of the temporal response, both the total 
population n t (a-c) and the spin redistribution 5nt (d), under 
resonant (solid curves) and non-resonant, with Alut p = 7r/3 
(dashed curves) and Alut p = 2n/3 (dotted curves), pumps. 
The excitation amplitudes, A, are indicated by arrows. 
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FIG. 3: Oscillations of Rabi flopping population in c-band, n, 
and the spin redistribution, 5n, versus excitation amplitude, 
A (the same curves as in Fig.l). 

citation case. The harmonic dependencies with periods 
A = 2tt and A = take place for the circular and 

linear polarized pumps respectively. By contrast, one can 
see a beating of n and Sn with A under the elliptically po- 
larized pump. Moreover, the spin orientation may change 
a sign for this case; once again, Sn = —n and Sn = for 
the circular and linear polarized cases respectively. 

Let us discuss the assumptions used in the above cal- 
culations. The main restriction is the consideration of 
the vicinity of the excitonic resonance without the con- 
tribution of two-photon biexcitonic resonance. Here we 
have applied the phenomenological initial conditions for 
Eq.(ll). This is valid if hAuj and Tl/t p are smaller than 
the exciton-biexciton splitting energy (around 3meV in 
InGaAs QDs) while a more detailed many-particle con- 



sideration is required for greater TlAlj and/or Tl/t p . We 
also have described the coherent response neglecting re- 
laxation because t p is shorter than the relaxation times. 
Next, the non-Markov equation (1) describes the inter- 
band generation of electron- hole pairs to second order ac- 
curacy, i.e. we have supposed that the energy \e\Ev$i_/ lu 
is much lower than the energy gap. Finally, the simple 
model of QD has been used in order to obtain Eqs.(8,12) 
and more complicated numerical calculations, e.g. sim- 
ilar to |l2j |. have to be performed in order to improve 
these simple formulas. 

In conclusion, we have considered the peculiarities of 
Rabi oscillations in a single quantum dot due to the spin- 
dependent ultrafast excitation by an elliptically polar- 
ized pulse. For the simple model of SQD with the spin- 
flip transitions forbidden we have found both the beating 
of Rabi oscillations under an elliptically polarized pump 
and the quenching of spin orientation due to spin-orbit 
splitting effect. A more detailed treatment of these ef- 
fects requires the consideration of complicated selection 
rules and to describe the exciton-biexciton resonance in 
the framework of a many-particle approach. What this 
paper does, however, is to determine the conditions for 
observation of the peculiarities mentioned. 
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